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Abstract 

Ferromagnetism with high Curie temperature T c , well above room temperature, and very small 
saturation moment has been reported in various carbon and boron systems. It is argued that the 
magnetization must be very inhomogeneous with only a small fraction of the sample ferromag- 
netically ordered. It is shown that a possible source of high T c within the ferromagnetic regions 
is itinerant electrons occupying a narrow impurity band. Correlation effects do not reduce the 
effective interaction which enters the Stoner criterion in the same way as in a bulk band. It is also 
shown how, in the impurity band case, spin wave excitations may not be effective in lowering T c 
below its value given by Stoner theory. These ideas are applied to CaBg and a thorough review of 
the experimental situation in this material is given. It is suggested that the intrinsic magnetism of 
the B2 and O2 dimers might be exploited in suitable structures containing these elements. 

PACS numbers: 75.10.Lp 75.50.Pp 
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I. INTRODUCTION 



There is great interest in the possible existence of ferromagnetism well above room tem- 
perature in materials which have no transition metal or rare earth metal components. Such 
high temperature^ ferromagnetism has been reported in various carbon systems such as 
fullerenes [l^ |2|j_ 3J, U, 0] and graphite J, 0, Q] and also in systems such as CaB 2 C2 



U\m c 



and CaB 6 |3, [l0( containing boron. The electrons involved in the ferromagnetism are s 
and p electrons, rather than d or /. The saturation moment in these systems is invariably 
very small and it is generally believed that the ferromagnetism is associated with defects or 
impurities. It is natural to suggest that the ferromagnetism may merely be due to magnetic 



impurities such as Fe but this appears to be ruled out experimentally in some cases 
It is also difficult to see how local magnetic moments associated with dilute Fe impurities 
could be coupled strongly enough to be ferromagnetically ordered above room temperature. 
In fact the same problem exists for local magnetic moments arising from any other form 
of defect. We believe that the most likely source of high temperature ferromagnetism is 
itinerant electrons occupying a narrow impurity band. 

In this paper we investigate how itinerant electron ferromagnetism in a narrow impurity 
band differs from the usual situation in the 3d band of transition metals. In the next section 
we discuss the criterion for ferromagnetism and the magnitude of the Curie temperature 
T c within Stoner theory. It is pointed out that the size of the effective on-site interaction 
parameter which appears in the theory is not limited by the width of the impurity band, 
as one might expect naively from Kanamori's T-matrix theory of electron correlation. 
A formal argument concerning this point is given in section 3. It is well-known that in 
transition metals, and in very weak itinerant ferromagnets such as ZrZn2, Stoner theory 
overestimates T c by a large factor. Instead T c is determined by low-lying spin fluctuations 
which are not included in Stoner theory . In section 4 we argue that this is not necessarily 
the case in impurity-band ferromagnets so that T c may in fact be close to its Stoner value. 

In all the carbon and boron-based ferromagnets the experimental situation is quite unclear 
with many conflicting results. To illustrate this we take the example of CaB 6 and in section 
5 many experiments on this system are reviewed. An attempt is made to build up a picture 
of the system based on the present theory of impurity-band ferromagnetism. In section 6 
we draw some conclusions. 
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II. STONER THEORY OF IMPURITY BAND FERROMAGNETISM 



A condition for ferromagnetism in an itinerant electron system is the familiar Stoner 
criterion 

I eff N(E F )>l, (1) 

where N {Ep) is the density of one-electron states per atom per spin at the Fermi level in the 
paramagnetic state and I e ff is an on-site interaction parameter. This criterion is not satisfied 
for most of the 3d metals and for none of the 4d metals. Even for the ferromagnetic metals 
Fe, Co and Ni it is only satisfied by a small margin (I e ffN (Ep) > 1.2) The reason 

for this is that the large on-site Coulomb interaction is reduced to an effective interaction 
I e ff ~ W/5, where W is the width of the <i-band, owing to a correlation effect [12]. Rather 
than experience the strong on-site interactions two electrons (or holes) avoid coming on the 
same site as far as possible and thereby increase their kinetic energy due to increased spatial 
confinement. Since N (Ep) ~ 5/W we have I e ffN {Ep) ~ 1. The situation is quite different 
when ferromagnetism occurs in a narrow impurity band. In this case N (Ep) ~ ni mp /Wi mp 
where n imp is the fraction of impurity atoms and the width Wi mp of the impurity band can 
be very small. Also, even if the bare interactions for sp electrons are strong enough for the 
Kanamori effect to operate, the kinetic energy increase due to increased spatial confinement 
of holes is governed by the full width of the valence band, not by the width of the impurity 
band. Thus 7 e // 3> Wi mp and this is shown formally in section 3 by considering the T-matrix 
for itinerant electrons in the impurity band. The Stoner criterion I e ffN (Ep) > 1 demands 
that n imp I e ff /W imp > 1 and this can be easily satisfied if W imp is sufficiently small. 
In Stoner theory the Curie temperature T c is given by 

Ieff J dE N {El) = 1 (2) 

where / (E) = {exp [(E — fi) /k B T c ] + l} -1 , /i is the chemical potential and N (E) is the 
density of states per atom per spin in the band. For simplicity we consider a half-filled 
impurity band with constant density of states N (E) = ni mp /Wi mp in the range — Wi mp /2 < 
E < Wi mp l 2. Here it is assumed that each impurity atom introduces one state of each spin 
into the band and provides one electron (or hole). Owing to the assumption of half-filling, 
(j, = at all temperatures so that Eq.(j2J gives 

k B T c = W imp / [4tanh -1 {W imp / I ef} n imp )\ . (3) 
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The right-hand side is a monotonic decreasing function of Wi 



so letting W%. 



imp 



we 



have 



ksT c < I e ffn imp / A. 



(4) 



To obtain a Curie temperature above room temperature, with a typical value of I e ff = 1 
eV, we require ni mp > 0.1. Assuming that there is complete spin alignment in the ground 
state, which for the rectangular band considered here occurs whenever the Stoner criterion 
n im P Ieff > W imp is satisfied, the saturation moment is therefore greater than 0.1 fi B / atom. If 
there is only partial spin alignment in the ground state the argument of section 4 supporting 
applicability of the Stoner model fails and low energy spin fluctuations are likely to reduce 
T c far below the Stoner value. We conclude that room temperature ferromagnetism and a 
uniformly distributed saturation moment much less than 0.1/iB per atom are incompatible. 

The observed magnetizations in graphite and Cgo are of the order 10~ 3 Am 2 /kg ~ 2 x 
10~ 6 hb per C atom but typically an order of magnitude larger in CaB 6 . We must conclude 
that the magnetization is very inhomogeneous with perhaps only a fraction 10~ 4 of the 
sample ferromagnetically ordered. Some evidence for this inhomogeneity in CaBg is discussed 
in section 5. 

III. THE EFFECT OF ELECTRON CORRELATIONS 

As we discussed in the previous section, renormalization of the effective Stoner interac- 
tion due to correlation effects is of crucial importance for understanding itinerant-electron 
ferromagnetism. Formally, this renormalization can be taken into account via the T-matrix 
approaaQQwMch is exact for the case of small electron (or hole) concentration but 
is qualitatively adequate for arbitrary band filling. We start with the general many-body 
Hamiltonian: 



H 



Ht + Hu 




(5) 
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where A = im are the site number (i) and orbital (m) quantum numbers, a =|, I is the 
spin projection, c + ,c are the Fermi creation and annihilation operators, H t is the hopping 
Hamiltonian, and the Coulomb matrix elements are defined in the standard way 



(12 H 34) = J drdrV4(r)^*(r> (r - r') ^ 3 (r)^ 4 (r'), (6) 



where we define for brevity Ai = 1 etc. Following Galitskii [15J] let us take into account the 
ladder (T-matrix) renormalization of the effective interaction: 



13 



T™' (in) 24^ = (13 \v\ 24)— i ( 13 \ v \ 57 > G 56 («*>) Gj's (in - iu) (68 T aa ' (in) 



ft 5I17S 



24 



(7) 

where uj = (2n + l)7TksT are the Matsubara frequencies for temperature fc^T = /3 _1 (n = 
0, ±1, ...). Using the spectral representation for the Green's function 

oo 
— oo 

substituting Eq.flHl) into Eq.(jZJ) and calculating the sum over Matsubara frequencies in the 



usual way 



161 ] one can rewrite Eq.(J7J) in terms of the spectral density matrix f)(x) : 



(13 |T 24) = (13H24)+^(13|u|57)(57|P(£;)|68)(68|T(S)|24), 

5678 

oo oo 

(57 |P (£7)1 68) = [ dx [ dy l ~ L {X) ~ f{y) Pb& (x) p 7S (y) (9) 
J J E-x-y 



— oo — oo 



where E is the real energy. In this section we neglect spin polarization since the effective 
exchange parameter I e ff in the Stoner criterion (jlj should be calculated in the paramagnetic 
phase. 

If we take into account only on-site Coulomb interaction (the Hubbard approximation) 
than the T-matrix turns out to be also diagonal in site indices and the matrix equation @ 
holds assuming that 1,2,... 8 label only orbital indices and f)(x) is local (on-site) spectral 
density. 

For simplicity, we will consider further the case with one orbital per atom and thus omit 
orbital indices. The effective interaction I e ff is given by T(£7) with E corresponding to 



the sum of energies of two occupied electron states [12j. Then from Q we have I e f 
v/(l — Pv) = —1/P for large v. Hence in the bulk-band case (Fig. 1(a)) I e ff ~ W 



NJE) 




FIG. 1: Schematic density of states for (a) a weak itinerant ferromagnet and (b) the ferromagnetic 
impurity band model 



The energy spectrum for the impurity band model is shown schematically in Fig. 1(b), with 
a broad main band (region I) and a narrow impurity band (region II). The figure shows an 
impurity band split from the bottom of the main band but the present considerations are 
identical for the case of holes in an impurity band split from the top of a valence band, 
as in the model of CaB 6 proposed in section 5A. To estimate different contributions to the 



6 



function P(E) one should take into account that 

dxd V _/_"V 1 ^2 
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/D (z) p (j/) ~ — In (^^— J ZbandZimp, (10) 



where 



■^imp — J dxp (x) , 
11 

J dxp(x) = 1 - Z imp (11) 



Zband 



1 

are total spectral weights of the impurity and main bands, respectively. 

To proceed further, one has to specify our model for the impurity band. Let us consider 
for simplicity the Slater-Koster model of single impurity, with non-degenerate band energy 
e(k), k being the quasimomentum, and on-site only impurity potential V (V < 0). Then 
the impurity site Green's function reads 

G m (E) = [F- 1 {E)-Vy 1 , 

*-w = *x:^7(kj (12) 

where N is the number of atoms in the crystal. The energy of the impurity localized state 
E and its spectral weight are determined by the equations 



VF (E t 



a j 



d 



dE 

For the case of shallow impurity levels 



e=e V*\F>(E ) 



(13) 



e = \E b - E \ < W, (14) 

where Ef, is the energy at the bottom of the main band, one can use the asymptotic forms 

1 f W \ 
F{E)^-\n[——) (15) 



W \Eh-E 
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and 



F(E) 



crit 



1 



+ B 



w 3 / 2 



(16) 



for two-dimensional (2D) and three-dimensional (3D) cases, respectively. Here V cr u is the 
critical potential for which an impurity state splits from the band and B is a dimensionless 
constant. We consider the 2D and 3D cases because the top of the valence band in CaBg 
has strong quasi-2D character, as discussed in section 5A. Substituting Eqs.([T5|). (JTfij) into 
Eq.ffT^j) one finds 



Taking into account Eqs. (jl3|) . (fT5|) and (jlfjjl one can eliminate V and rewrite this estimation 
in the final form 



taking BV 2 rit /2 ~ W 2 . One can see that for the case of shallow levels the spectral weight 
of the impurity state is small. This estimation holds, at least in order of magnitude, also 
for many-impurity case when the impurity band is formed. Thus one can conclude that all 
contributions (fTUj) to -P (E) are much smaller than l/Wi mp so that J e // 3> Wi mp . Thus the 
Kanamori renormalization of the effective Stoner parameter is not limited by Wi mp as it is 
limited by W in the bulk case. . 

IV. SPIN WAVES IN THE IMPURITY BAND MODEL 

In Fig. 1 we contrast two types of itinerant electron ferromagnet with small magnetization 
in the ground state. Case (a) shows schematically the situation in a very weak itinerant 
electron ferromagnet such as ZrZn2. The Stoner criterion is only just satisfied and a small 
exchange splitting between the majority (|) and minority (|) spin bands leads to the Fermi 
level Ep being positioned slightly differently in the two bands. Consequently, the majority 
spin Fermi surface is slightly larger than the minority spin one and the volume of k space 
between the two Fermi surfaces is proportional to the ground state magnetization. The Curie 
temperature is determined by collective spin fluctuations rather than the single particle 
excitations of Stoner theory [13]. Well-defined spin waves play a negligible role because 
they are confined to a very small region of the Brillouin zone. However for sufficiently small 




(17) 




(18) 
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wave-vector q they exist with energy Dq 2 . The spin wave stiffness constant D is small, being 
proportional to the magnetization 18[ . 

The situation is different in case (b) where the small magnetization is associated with 
complete spin alignment of a low density of carriers in an impurity band. The dominant low- 
lying spin-flip excitations are now long-wavelength spin waves and we proceed to estimate 
D in this case. It is shown below that as long as complete spin alignment is maintained in 
the ground state D no longer tends to zero with the magnetization. 

For simplicity we consider a one-band model with the disorder leading to the impurity 
band treated within the coherent potential approximation (CPA). We suppose that near 
the bottom of the band the Bloch-state energy e (k) in the absence of disorder and local 
exchange splitting is given by 

h 2 k 2 

£ (k)=a + — . (19) 

Within the CPA the effect of disorder and exchange is contained in uniform potentials, for 
an effective medium, which depend on spin and energy. These potentials, or self-energies, 
So- (E) for spin a, would have to be determined self-consistently for a given detailed model, 
but here it is sufficient to assume that they give rise to exchange-split impurity bands as in 
Fig. 1(b). 

The CPA Green's function for spin a is of the form 

Gl (E) = [E- e (k)-^(E)]- 1 (20) 



with spectral representation (cf Eq.( 



kV ; J E-x + irj J E-x + iT] ' y J 

1] — > +0. Here the spectral functions pi mp and pband refer to the impurity band and main 
band respectively. 

The spin wave stiffness constant is given by 

E F 

D = i ^ f dE^ \G{ (E) - Gi (E)} 2 |V k e (k) | 2 (22) 



where n a is the number of carriers of spin a 



19, 201. This result was first derived for 



ferromagnetic alloys within the random phase approximation but is valid within any local 



approximation, where the self-energy is a function of energy only 



21 



In the case of complete spin alignment considered here, = and (E) is real for 
E < Ep. By using Green's theorem on the k sum we find 

D = J-9 f dEj2\-Gl(E)Vle(k)-2Gi(E)Gi(E)\V k e(k)\ 2 

-oo k 
E F 

= ^ / dEj2[pi(E)Ke^) + 2 P l(E)Gi(E)\V k e(k)\ 2 ]. (23) 

-oo k 

Noting that Vj^e (k) /6 = h 2 /2m*, and using the spectral representation of G^., we may write 

D = D + D 1 (24) 



where 



D = h 2 /2m* (25) 



Dl = _L| d£^pt(£) | |V k e(k)| 2 . (26) 



— oo 



To evaluate the energy integral in .Di we approximate the denominator E — x in Eq. (j2T?j) by 
Dot — x > an d further approximate it by —A = E ^ — E ^ in the integral involving p^. imp (x) . 
Here E 0a is the energy at the centre of the a spin impurity band and A = Eq± — E^ is the 
splitting between the two impurity bands. Thus D\ becomes 



where (n£) = J dEp%_ (E) . To determine the k dependence of ( n^) and p^ imp it is sufficient 

— oo 

to consider a single impurity with the on-site potential V, as in the previous section. Then 

Gk (E) = - 1 - , + — 2 (28) 
1 ^ E-e(k) [£-e(k)f 1 ' 

with t (E) = V/ [1 — V^F (£")] (cf. Eq. (jl2|) \ For the case of a single impurity the contribution 
of the impurity level to the spectral density reads, similar to Eq. ffTB*]) : 



'^ Gk{E) -[E -e(k)fF>(E ) 



3G k (E) = - 2 5 (E - Eo) . (29) 
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For a small density of impurities the ^-function broadens to the shape of the impurity band 
Pi mp (E) and we may write 

■ (30) 

Ep oo 

Now, n-f = / dEp\ mp (E) and N imp = J dEp% mv [E] , where N imp is the total number of 

— oo — oo 

impurities and it is assumed that each impurity contributes one state to the impurity band. 
Using these equations, together with Eq. (|HU|) . and assuming for simplicity that Eq^ ~ Eqi ~ 
Eq (i.e. A^Bj- Eq), we obtain the contribution to D\ of the first term in curly brackets 
in Eq. (|27p in the form 

3m- A \^ [E - e (k)] 2 j ^[£ - e (k)] 4 1 ' 

To estimate the k sums in this expression we assume a constant density of states C per 
atom in the band e (k) . The first term in Di, given by 1)31 J), the becomes 

^ 2 AW ^ 2 ft 2 

2m*9/ e// n T C 9 2m* 1 J 

if n T ~ A^ mp and I eff C ~ 1. 

To estimate the second term in D\ we assume that the bulk of the band is only weakly 
perturbed by the disorder so that p^band (E) ~ 8 (E — e (k)) . Combining the result with 
Eq.fjHU), and using Eq. (|2*3|) . we obtain 

D = J D + J Di^^ = ^ . (33) 
9 2m* 9 

This is a crude estimate for an oversimplified model. The main point is that, as long as 
complete spin alignment is maintained, D is independent of the ground state magnetization. 
This contrasts with the standard case of very weak itinerant electron magnetism where D 
is proportional to the magnetization Q|. 

A realistic estimate of Do, the first term in D, may be obtained within a multi-orbital 
tight-binding model (e.g., Ref. 3|). To make the calculation specific we consider the case of 
CaBg where, as discussed in the next section, we propose that ferromagnetism arises from 
the complete spin alignment of holes in an impurity band just above the valence band. The 
valence band is assumed to be formed largely from boron 2p orbitals and the Hamiltonian 



11 



is of the form 



H = H + Hi , 

H ° = EE^'(k)4 CT v kff , (34) 

k flfj,'cr 

where c^ ko . creates an electron of spin a in a Bloch state of wave-vector k formed from 
orbital /i. Here Hi contains on-site Coulomb and exchange interactions and an additional 
one-electron term representing spin-independent diagonal disorder. For a cubic crystal, such 
as CaB 6 , 

D » = 6 fr_ n) E E V ^ (k) (4 ffVkCT ) (35) 

where the expectation value is evaluated in the ferromagnetic ground state. If we assume 
only nearest-neighbour hopping, and the on-site orbital energy is taken as zero, V k V^ M ' (k) = 
—R 2 V flfl / (k) where R is the nearest-neighbour distance. (We neglect the slight difference in 
distance between boron atoms within the unit cell of CaB 6 and between unit cells). Hence, 
from Eq.(j2SJ), 

A, = ^ S (Ho) = ^R 2 W, (36) 

6 (n T - n;) 12 

where W is the width of the valence band. The last equality follows since (H ) = for a 
completely full band and rif — n; holes are removed from an impurity band formed from 
states at the top of the band with energy W/2. D is an upper bound on D since the second 
term in D is ne gati ve. For CaBg the B-B distance R — 1.7 Aand the width of the valence 
band W = 9 eV |23|]. Hence D = 2167 meV- A 2 =3.47xl(T 39 Jm 2 . At low magnon density 
the magnons may be regarded as non- interacting bosons with energy Dq 2 . The number of 
magnons excited per unit volume at temperature T is then 0.0586 (kBT/D) 3 ^ 2 . The volume 
per boron atom in CaB 6 is 11.88 A 3 so that at T = 600 K with D = D the number of 
magnons excited per boron atom is 2.57 x 10~ 3 . The corresponding reduction in moment 
per boron atom is 5.14 x 10~ 3 /i£ which is much less than the minimum T = saturation 
moment of 0.1/i£ discussed in section 2. Thus if D = Dq spin waves are not the dominant 
excitations controlling T c and the Stoner estimate of this quantity should be realistic. If, 
however, D is strongly reduced from D , to D /5 say, the number of magnons excited is an 
order of magnitude larger and T c could be reduced considerably from its Stoner value. Thus 
to obtain a Curie temperature well above room temperature we require D > -Do/2, say. 
Whether this is possible can only be decided by realistic calculations of D in a multiorbital 
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impurity band model. This can certainly be done, but only if one knows the nature of 
the impurities or defects producing the impurity band. No reliance can be placed on the 
calculation for a simple one-band model discussed earlier in this section which gave the 
estimate D ~ D /9. 

To conclude this section we point out why it should be much easier to achieve high tem- 
perature ferromagnetism in the impurity band system than in diluted magnetic semicon- 
ductors such as (Ga,Mn)As. In our impurity band model itinerant electron ferromagnetism 
is achieved through the high density of states in a narrow impurity band, which enables 
the Stoner criterion to be satisfied. In the diluted semiconductor ferromagnetic order is 
obtained by an exchange interaction between the local Mn moments which is mediated by 
spin polarized holes in the semiconductor valence band. In this case Eq. (J35j) for D Q is still 
valid but the factor — includes the moment of the Mn atoms. Thus, for the same 
number of spin aligned holes in the band (or impurity band) in both models, and similar 
bandwidth W and interatomic distance R, the relative value of D$ is controlled by the 

— rii factor. Hence in the diluted ferromagnetic semiconductor Do is reduced by a factor 
Mband/ (Mb ani i + M/ oc ) compared with the impurity band case. Here Mb an d and M/ oc are the 
contributions to the saturation magnetization of band carriers and localized moments re- 
spectively. With Mi oc = 4{xb for Mn atom this factor can be small. The role of spin waves in 
reducing the Curie temperature of diluted ferromagnetic semiconductors has been discussed 
by Jungwirth et al [24]. The reason why local moments reduce D so strongly is that they 
provide no exchange stiffness, this being entirely due to the band carriers. It is like hanging 
weights on a spring without increasing its stiffness. The oscillation frequency is reduced. 



V. APPLICATION TO CAB f 



In this section we review most of the experimental data on CaBg, doped and undoped, 
and try to build up a theoretical picture based on the impurity band model. Most other 
theories are based on homogeneous models where doping merely changes the number of 
electrons in the system. The first such model ^, [2^, [2(^ assumes the system behaves as a 
low-density electron gas which is predicted to become ferromagnetic for density less than 
about 2 x 10 20 cm~ 3 (r s > 20). Even favourable estimates of T r 291 are far below room 
temperature. Another early model 

0, H H Q 

is that of a doped excitonic insulator 
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which was predicted long ago to be ferromagnetic |3jj. Problems associated with lack of 
nesting between electron and hole Fermi surfaces and with structural instability have been 
discussed j^J. Hotta et al ,32] show that the excitonic state is unlikely to be stable and also 
conclude that ferromagnetism in CaB 6 may require the local removal of cubic symmetry by 
defects. 

Young et al first observed high-temperature ferromagnetism with small saturation 
moment in La-doped CaBg, SrBg and BaBg. In Cai-^La^Be with x = 0.005 the saturation 
moment M s — 3.5 x 10~ 4 /is/unit cell and the Curie temperature T c = 600K. Subsequently 
Ott et al [loj found M s = 2 x lO^s/unit cell and T c = 900K for x = 0.01. According 
to the discussion of section 2 the magnetization must be very inhomogeneous with only a 
fraction 10~ 4 — 10~ 3 of the sample ferromagnetically ordered. 

Young et al jft] found that ferromagnetism does not appear for x > 0.02. This has been 



confirmed recently by Cho et al [ll| in Cai_ x La x B 6 crystals with x = 0.03 and 0.04 grown 
using boron of 99.9% purity (3N). They did find ferromagnetism in nominally stoichiometric 
CaBe, as well as in Cai-ajLa^Be with x = 0.005, 0.01 and 0.02. All these ferromagnetic 
crystals were grown using 3N boron. Cho et al |Tjj found no magnetism in boron-deficient 
CaB 6 (3N) and in no single crystal of CaB 6 (stoichiometric, boron rich or boron deficient, 
or La-doped) made using boron of 99.9999% purity (6N). Furthermore the temperature- 
dependent resistivity p(T) of CaB 6 (6N) exhibits semiconducting behaviour, corresponding 
to an energy gap much smaller than the 168 meV deduced by Vonlanthen et al |23J] from 
their p(T) data, whereas CaBg (3N) is metallic with p(0) a factor 10~ 4 smaller than in CaBg 
(6N). Hall effect measurements on ferromagnetic undoped CaBg crystals invariably indicate 



that electrons are the dominant carriers 



34j . From the work of Cho et al 11] it appears 



that ferromagnetism is associated with charge carriers introduced by impurities in 3N boron. 
The removal of ferromagnetism by overdoping with La, which adds electrons to the system, 
suggests that ferromagnetism is associated with holes in the valence band, or in an impurity 
band close to the valence band. In a ferromagnetic region, with 0.3 aligned hole spins per 
unit cell, say, 2% La doping would remove 0.02 holes per unit cell which might be enough to 
move the Fermi level out of a sharp peak in the density of states so that the Stoner criterion 
is no longer satisfied. Clearly if the magnetization were uniform the critical concentration 
for uniform La doping in this picture would be more like 0.002%. Boron deficiency would 
also be expected to remove holes from the valence band. However Hall effect measurements 



14 



conduction band 



impurity 
bands ♦ 



valence band 




N(E) 



FIG. 2: Sketch of the proposed electronic structure of ferromagnetic C&Bq. 

34I ]. particularly those 35] on a ferromagnetic sample with p(T) very similar to the CaB@ 
(3N) crystal, indicate that electrons are present in the conduction band. This suggests 
that impurities or defects arising from 3N boron lead to a small band overlap. The major 
impurities in 3N boron were identified as C and Si, with negligible amounts of magnetic 
impurities such as Fe jll|. 



A. An impurity band model of ferromagnetism in CaBg 

The first assumption of the model, which seems essential for reconciling high T c with 
very small average moment per atom, is the inhomogeneous nature of the magnetization. 
Ferromagnetic regions are assumed to have a moment of order 0.1 \xb per B atom, about 
1000 times the average value. As discussed in the previous section this is consistent with 



the high La doping concentration required to suppress ferromagnetism. Terashima et al £ 
invoke spatially-inhomogeneous magnetization in Cao.ggsLao.oosBe to try to understand their 
susceptibility data. In ferromagnetic resonance measurements on the same system Kunii 
3?| finds a magnetization which is three orders of magnitude larger than the value observed 
in static measurements. This could be consistent with our picture. 

The second assumption of the model is that the saturation moment within a ferromagnetic 
region arises from complete spin alignment of holes in a narrow impurity band just above the 
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valence band. Electrons are also present in an overlapping conduction band whose states, to 
a good approximation, do not mix with those of the valence band. The situation is shown in 
Fig. 2. Transport is dominated by the conduction band which may not be much influenced 
by the magnetic inhomogeneity of the valence band. 

To proceed further it is necessary to describe the nature of the CaB 6 system. The crystal 
is cubic, with Ca 2+ ions forming a simple cubic lattice. Each cubic unit cell is occupied by 
an octahedron of 6 boron atoms, but the distance between neighbouring B atoms in adjacent 
cells (~ 1.67 A) is considerably less than the distance between B atoms in one octahedron 
(~ 1.76 A) [23]. Thus the structure could be described as a Ca simple cubic lattice with B 2 
dimers positioned at the centre of each square face and oriented perpendicular to that face. 
Dimerization is a marked effect, the interatomic spacing in the isolated B 2 dimer being 1.59 
A, and it plays a significant role in our theory, as will be seen later. Clearly the interatomic 
spacing 2ua for dimers in the lattice can be varied uniformly without changing the cubic 
symmetry or the lattice constant a. Calculations by Massidda et al j3] predict that the 
band gap in CaB 6 reduces by Q.heV for a 3.5% decrease in u with constant a. The bottom 
of the conduction band, at the point X in the zone, falls while the top of the valence band 
at X rises. Negative band-gap corresponds to band overlap. A change from a very small 
band gap in the CaB 6 (6N) crystals of Cho et al [ll| to band overlap in CaB 6 (3N) might 
occur by a slight contraction of the dimers induced by impurities in 3N boron. This impurity 
effect would not lead to uniform contraction and an impurity band might form just above 
the valence band. This contraction effect might be produced by many different impurities, 
so the mechanism leading to ferromagnetism is not confined to the CaBg (3N) crystals of 



Cho et al [ll|. 

The formation of a shallow impurity band is favoured by the quasi-two-dimensional nature 
of the band-structure near the top of the valence band at the X point. All existing band 
calculations find a very flat band along XT which leads to an almost discontinuous rise in 
the density of states at the top of the valence band. However band calculations discussed 
in section 5B differ widely in their prediction of the sign and magnitude of the direct band 
gap at the X point in CaB 6 . 

The impurity band picture has similarities with the situation in ferromagnetic transition 
metals such as Ni and Co. The role of the transition metal d band is played by a narrow 
impurity band formed from boron p orbitals and in both cases there is a broad overlapping 
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conduction band, this being formed largely from calcium d orbitals in the CaBg case. Mag- 
netism in the transition metal d band is not unexpected since in the isolated atom the spins 
of electrons in degenerate d orbitals are aligned by the Hund's rule mechanism. An analo- 
gous situation in the boron p band becomes clear when one remembers, as discussed above, 
that the boron lattice may be regarded as formed from B 2 dimers, and that the isolated B 2 



38, 



39|. The two p electrons in B 2 occupy 



dimer has a ground state with total spin S = 1 
degenerate ir x , n y bonding orbitals, the z axis being the axis of symmetry of the dimer, and 
their spins are aligned according to Hund's rule. The only other magnetic second-row dimer 
is 2 , the relevant orbitals now being the corresponding antibonding orbitals. 



B. Further discussion of theory and experiment 

The picture described above is controversial for several reasons. The first controversial 
point to consider is whether CaBg is a semimetal (semiconductor), with a very small negative 
(positive) band gap, or whether it is a semiconductor with a gap of order 1 eV. Following 
the pioneering band-structure calculation of Hasegawa and Yanase (^(J, using the local 
density approximation (LDA), Massidda et al found that within LDA the dimerization 
effect discussed above led to a small band overlap. However the difference between small 
negative and positive band gap is within calculational error. The top of the valence band 
near X is formed from B bonding p orbitals and the bottom of the conduction band at 
X is dominated by Ca d orbitals. These states have a different parity on symmetry lines 
and the bands therefore cross each other there. Although mixing occurs at general points 
in the zone these two bands may therefore be considered approximately as independent of 
each other, even when they overlap. Rodriguez et al |4l] have made similar calculations for 
SrBg and find quite detailed agreement with optical and transport data |42| which support 
the semimetallic picture. Further strong support for this picture comes from more recent 
data on reflectivity and Hall resistivity. To obtain a reflectivity which relates to the 
bulk band-structure it was necessary to file the surface of the as^rown crystal to a depth of 



lOfim and then polish it. According to our discussion of Ref.|llj. CaB 6 and SrB 6 crystals 
are expected to be semimetallic if the boron used to make them is of substantially less 
than 6N purity. Hall et al |4jJ have made de Haas - van Alphen (dHvA) measurements on 
Cai-aXazBg, with x = 0,0.0025,0.005 and 0.01, and on SrBg. The La-doped samples were 
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ferromagnetic, except for x = 0.01 where in any case no dHvA oscillations were seen owing 
to impurity scattering. Orbits attributed to electron and hole pockets were seen in CaB 6 and 
SrB 6 but the hole pockets dropped out as electrons were added by doping with La. Electron 
and hole pockets have also been seen in the related divalent hexaboride EuB 6 by dHvA 4^ 
and Shubnikov - de Haas experiments The calculations and data so far discussed in this 
paragraph are reasonably consistent with the pictured we proposed in section 5A. However 
Denlinger et al |47[ report angle resolved photoemission spectroscopy (ARPES) data, on 
cleaved (001) surfaces, showing an A-point band gap of about 1 eV for SrB 6 and EuB 6 . 
The Fermi level is located just above the bottom of the conduction band. These authors 
interpreted the apparent discrepancy between their ARPES data and the bulk semimetal 
picture in terms of an electron-rich surface region, maybe 60 Athick, within which the Ap- 
point band gap varied rapidly due to variation of the dimerization parameter u between 
bulk and surface. They postulated that broken dimer bonds at the cleaved (001) surface 
would allow octahedra near the surface to contract, thus increasing u and consequently the 



band gap. Subsequently Denlinger et al [48, 49] reinterpreted their large observed gap as 
characteristic of the bulk. This was largely in response to a band calculation for CaBg by 
Tromp et al j3] who used the GW approximation to obtain a bandgap of 0.8 ± O.leV. 
The GW approximation generally yields larger, and more reliable, bandgaps than LDA. 
Denlinger et al support their new conclusion with results of X-ray absorption and emission 
spectroscopy although the band edges in these data are not clear-cut. Very similar ARPES 
data are reported by Souma et al |35J]. A more recent band calculation for CaB^ using 

nn 

the GW approximation is that of Kino et al |51l . I52| . They find that over most of the 
zone the gap between conduction and valence bands is enlarged compared with the LDA 
case, as in Ref. [50], but that a small band overlap nevertheless occurs at the X point. 
The subtle reason for this exceptional situation is analysed and it is suggested that the 
discrepancy with Ref. 0] may be due to an improper treatment of the core contributions in 
the pseudopotential plane-wave method by Tromp et al. A similar pseudopotential method, 
but with a non-local weighted density approximation (WDA), was recently used by Wu et 
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point out that the 



al |53( to find the same band gap as that of Ref. [50j . Kino et al 
difference between a very small negative and positive band gap is within their calculational 
error. The question of semimetal, or small band gap semiconductor, versus leV band gap 
semiconductor in the bulk is clearly not completely settled. However there is strong evidence 
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for the former scenario 143 



44L l5ll l52| and we have adopted this viewpoint here. 
The second controversial point we must consider is whether the small-moment ferromag- 
netism in these systems arises from spin alignment of charge carriers within ferromagnetic 
regions of the bulk or is associated entirely with magnetic impurities near the surface. It is 
also quite possible for these two sources of magnetism to coexist. The standard method of 
growing CaBg crystals uses an Al flux and Otani and Mori (54], |2| attribute ferromagnetism 
to Fe on the surface, the Fe being introduced as an impurity in the Al flux. They found 
that ferromagnetism disappeared completely from crystals of CaB 6 and Cai_ x La x B 6 after 
being kept in HC1 solution for a few hours. They conclude that Fe was removed from the 
surface during the HC1 treatment and that no bulk ferromagnetism exists in these samples. 
However many workers have used the same Al flux to grow a number of crystals, some ferro- 
magnetic and some not. For example ferromagnetism never exists in overdoped Cai_ x La x B6 



with x > 0.03 



11] . Cho et al ll| found that no crystal made with 6N boron was ferro- 



magnetic, even though the same Al flux was used as for ferromagnetic 3N boron samples. 
However in some of the latter samples it was found that chemical etching partially removed 
the magnetism, although no details are given. Thus bulk ferromagnetism and surface im- 
purity ferromagnetism may coexist in these samples. Bennett et al [34j also find evidence 
for such coexistence in a ferromagnetic sample of CaBg. The saturated moment M s was 
reduced by 47% after a surface layer of thickness 6000 Awas removed by etching, but no 
further reduction occurred on removing an additional 6000 A. These workers also found that 
in 16 different undoped CaBg crystals there was no correlation between M s and the electron 
density deduced from Hall effect measurements. This is consistent with our picture that the 
bulk part of M s is determined by the number of holes in the valence band, or an associ- 
ated impurity band, whereas the Hall effect is dominated by conduction electrons. Young 
et al [56j report measurements of M s for CaB 6 grown from Al flux which was deliberately 
contaminated with a wide range of concentrations of added Fe. There was no dependence 
of M s on Fe concentration, from which the authors conclude that alien Fe-B phases are 
probably not the source of ferromagnetism in aluminium-flux-grown single crystals. This 
contrasts with the behaviour of CaB 6 sintered powders for which Matsubayashi et al 57 1 
offer convincing evidence that the observed ferromagnetism is entirely due to FeB and Fe?B 
phases on the surface, the Fe originating in the crucible used for synthesis. Young et al 56| 
propose that in aluminium-flux-grown crystals Fe impurities incorporated in the bulk are 
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responsible for high-temperature ferromagnetism. It is assumed that the concentration of Fe 
is about 0.1 atomic per cent, being determined by off-stoichiometry of the C&B$. It seems 
that such a low concentration of Fe, considered purely as a magnetic defect, could only lead 
to ferromagnetism with a very low T c . In fact Young et al report that inclusion of Co or 
Ni in the flux, in addition to Fe, leads to small-moment ferromagnetism with T c < 10K. 
They propose that Co and Ni enter CaBg in preference to Fe but do not explain why Fe 
should lead to a T c 100 times larger than is produced by Co or Ni. It is possible that Fe is 
one of the defects which lead to a partially occupied impurity band just above the valence 
band, whereas Co and Ni do not. We have argued that such an impurity band can lead to 
high-temperature ferromagnetism; this is quite independent of whether or not the defects 
responsible for the impurity band are magnetic impurities. Meegoda et al j3] argue against 
any significant concentrations of Fe in the bulk on the basis of depth profiling with Auger 
electron spectroscopy; for a CaB 6 crystal grown in a Al flux they find that Fe impurities are 
confined to a surface region with a depth of a few microns. 

From the above discussion we conclude that high-temperature ferromagnetism of CaE>6 
and Cai_ x La x B6 crystals exists as an inhomogeneous bulk effect. However in as-grown 
crystals there is frequently a contribution from Fe impurities in a surface region which can 
be eliminated by suitable surface treatment. The bulk ferromagnetism is clearly associated 
with defects or impurities but, as in the case of impurities arising from use of 3N boron in 
Ref . [ll| , they need not be intrinsically magnetic. Jarlborg 59, 60] looked for ferromagnetism 
associated with various defects in StBq by means of LDA band calculations. The point 
defects were centred in 2x2x2 or 3x3x3 periodically continued supercells, corresponding 
to defect concentrations x = 0.125 and 0.037 respectively. The point defects considered 
were La, In and Al impurities, and a vacancy, replacing one Sr atom. For the La case 
a moment of order 0.1/i£ per La impurity was found for both sizes of supercell. However 
similar calculations for La in CaBg by Monnier and Delley [6J| , using the generalized gradient 
approximation (GGA), a variant of LDA, and a 3 x 3 x 3 supercell, resulted in zero moment. 
Jarlborg emphasizes that for ferromagnetism it is essential that the Fermi level falls within 
an impurity band, which in the case of La is associated with the conduction band. No 
magnetism is found in the In, Al or vacancy case. For In and Al doping with x = 0.125 
Jarlborg finds an almost rigid upward shift of the energy bands relative to the Fermi level, 
by about O.SeV^ in the case of In, so that holes occupy the valence band. A similar effect also 



20 



occurs for the Sr vacancy. For lower doping levels presumably no rigid shift would occur and 
the holes might occupy a narrow impurity band just above the valence band. Jarlborg does 
not consider any lattice relaxation around the impurity. Monnier and Delley [61] find a large 
magnetic moment of 2.36yU# associated with removal of a complete boron octahedron; this 
moment is reduced to 1.32/xb if the B 6 octahedron is replaced by Ca and is lost completely 
if E>6 is replaced by La or Al. They propose that the large moments associated with B$ 
vacancies, created during growth due to kinetic effects, are the source of ferromagnetism. 
Fisk et al (3] take up this suggestion and estimate the density of such vacancies to be 
10~ 4 — 1CT 3 per unit cell. However no explanation is given of how such a low concentration 
of local moments could have a Curie temperature in the range 600 — 900K (see also Ref. 63j). 

From the above discussion of Jarlborg's works it appears that Al substituting for Ca or 
Sr, or a Ca or Sr vacancy, are possible sources of the impurity band required for ferromag- 



netism. Terashima et al |36( deduce from a low temperature resistivity anomaly, also seen 
by Vonlanthen et al fj^j], that their CaBg samples were contaminated with Al granules. It 
therefore appears that Al does enter the bulk during Al flux growth and might also sub- 
stitute for Ca. Hall et al [44] find no evidence of Al inclusions in their samples but Al as 
a substitutional impurity is presumably not ruled out. However Monnier and Delley 61 1 
find that the formation energy of a Ca vacancy and an Al substitutional impurity are both 
about 5eV, so that the defect concentrations would be negligible if growth took place under 
conditions of thermal equilibrium. This is probably not the case since it is generally be- 
ieved that hexaborides, including CaB 6 , have a tendency to self-doping by metal vacancies 

[t might be that certain non-magnetic impurities, including those 
in 3N boron which are apparently essential for ferromagnetism according to Cho et al 
facilitate the formation of Ca vacancies or Al substitutional impurities by kinetic effects 
during growth, as suggested by Monnier and Delley [6l| in connection with B 6 vacancies. 
On the other hand, we have already described how any impurity leading to more pronounced 
boron dimerization locally could lead to the desired impurity band. The effect of Fe impu- 
rities in the bulk on the electronic structure is unknown. The same can be said of C and Si 
impurities which are probably introduced using 3N boron in the work by Cho et al [ill ]. 

We suggest that there are many ways in which a narrow impurity band can be formed 
just above the valence band and we believe that this is essential for ferromagnetism. The 
extremely inhomogeneous nature of the ferromagnetism which we propose might arise in 
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two ways. One possibility is that the impurities or defects responsible for the impurity band 
might be distributed very inhomogeneously. A second possibility is that the impurity band 
is widespread but the Stoner criterion for ferromagnetism is only satisfied in regions with a 
particularly favourable impurity configuration. 

VI. OUTLOOK 

A ferromagnetic semiconductor with a Curie temperature T c well above room temperature 
would have great potential for use in spintronic devices. A well-explored route towards this 
objective involves Mn doped III-V compounds but the highest T c obtained is well below 
room temperature. The reason for this failure to achieve high T c is discussed at the end of 
section 4. The essential point is that although ferromagnetic order is induced by the large 
Mn moments they do not contribute to the exchange stiffness, which is provided by the spin- 
polarized band carriers. In fact the Mn moments reduce the spin wave stiffness constant by a 
large factor. To achieve high T c it seems essential that the band carriers themselves produce 
ferromagnetic order by their mutual interaction, without the aid of local magnetic moments. 
This requires a high density of states at the Fermi level, to satisfy the Stoner criterion, and 
this is most likely to be achieved in a narrow impurity band. However the estimate of T c in 
section 2 indicates that room temperature ferromagnetism is unlikely unless the number of 
completely spin-polarized carriers in the impurity band exceeds 0.1 per bulk atom. It is not 
clear whether this can be compatible with a sufficiently narrow impurity band. Assuming 
that this scenario is possible, at least in some regions of the sample, it is argued in section 
4 that spin wave excitation may not significantly reduce T c from its Stoner value. If this 
mechanism is responsible for the ferromagnetism observed in some carbon systems and in 
CaB 6 , the very small observed magnetization indicates that ferromagnetic order exists only 
in a very small fraction of the total volume of the crystal. Clearly, before any spintronic 
application can be envisaged, it is necessary to find a way of producing the conditions for 
impurity band ferromagnetism more uniformly throughout the specimen. 

In section 5A we pointed out that the CaB 6 structure may be regarded as built from B 2 
dimers. The isolated dimer is magnetic, with a moment of 2 \xb ', of course this large moment 
does not survive in the CaB 6 structure where the bonding levels broaden to form the valence 
band as electrons hop between bonding orbitals of adjacent dimers. It would be interesting 
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to find a boron structure with more pronounced dimerization, and narrower bands formed 
from the bonding orbitals, than in CaB 6 . If such a structure exists its magnetic properties, 
doped and undoped, could be interesting. The same can be said of any structure containing 
O2 dimers, linked by hopping between the antibonding orbitals of the dimer in this case. 
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